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ABSTRACT 



We describe a systematic program aimed at identifying and characterizing candidate high-mass proto- 
stellar objects (HMPOs). Our candidate sample consists of 69 objects selected by criteria based on 
those established by Ramesh & Sridharan (1997) using far-infrared, radio-continuum and molecular line 
data. Infrared- Astronomical-Satellite (IRAS) and Midcourse-Space-Experiment (MSX) data were used 
to study the larger scale environments of the candidate sources and to determine their total luminosities 
and dust temperatures. 

To derive the physical and chemical properties of our target regions, we observed continuum and 
spectral line radiation at millimeter and radio wavelengths. We imaged the free-free and dust continuum 
emission at wavelengths of 3.6 cm and 1.2 mm, respectively, searched for H2O and CH3OH maser emission 
and observed the CO J = 2 — > 1 and several NH3 lines toward all sources in our sample. Other molecular 
tracers were observed in a subsample. 

While dust continuum emission was detected in all sources, most of them show only weak or no emission 
at 3.6 cm. Where detected, the cm emission is frequently found to be offset from the mm emission, 
indicating that the free-free and dust emissions arise from different subsources possibly belonging to the 
same (proto)cluster. A comparison of the luminosities derived from the cm emission with bolometric 
luminosities calculated from the IRAS far-infrared fluxes shows that the cm emission very likely traces 
the most massive source, whereas the whole cluster contributes to the far-infrared luminosity. Estimates 
of the accretion luminosity indicate that a significant fraction of the bolometric luminosity is still due to 
accretion processes. The earliest stages of HMPO evolution we seek to identify are represented by dust 
cores without radio emission. 

Line wings due to outflow activity are nearly omnipresent in the CO observations, and the molecular 
line data indicate the presence of hot cores for several sources, where the abundances of various molecular 
species are elevated due to evaporation of icy grain mantles. Kinetic gas temperatures of 40 sources are 
derived from NH3 (1,1) and (2,2) data, and we compare the results with the dust temperatures obtained 
from the IRAS data. 

Comparing the amount of dust, and hence the gas, associated with the HMPOs and with ultracompact 
Hn regions (UCHns) we find that the two types of sources are clearly separated in mass-luminosity 
diagrams: for the same dust masses the UCHn regions have higher bolometric luminosities than HMPOs. 
We suggest that this is an evolutionary trend with the HMPOs being younger and reprocessing less 
(stellar) radiation in the IR than the more evolved UCHns regions. 

These results indicate that a substantial fraction of our sample harbors HMPOs in a pre-UCHu region 
phase, the earliest known stage in the high-mass star formation process. 

Subject headings: stars:formation - stars: massive - ISM: hot-cores 



High-mass stars exert a decisive influence on the appear- 
ance and evolution of galaxies. Throughout their life cy- 
cle, they inject significant amounts of energy and momen- 
tum into their environments through stellar winds, molec- 
ular outflows, ultraviolet (UV) radiation and eventually 
supernova explosions. They also influence the formation 
of nearby low-mass stars, most of which are formed in 
dense clusters heavily affected by a few massive stars, the 
Trapezium cluster in Orion being a well studied nearby 



1. INTRODUCTION 



example (Hillenbrand & Hartmann 1998). Clearly, ques- 
tions concerning the number, birth-rates, distribution, and 
evolutionary timescales of massive stars are of critical im- 
portance for building a general picture of Galactic evolu- 
tion, in addition to their being interesting in themselves. 
While a scenario exists for the formation of low mass stars, 
much less is known about the processes involved in the for- 
mation of high-mass stars, which we take to be stars with 
M> 8Mq (corresponding to stars of spectral type earlier 
than B2) not having an optically visible pre-main-sequence 
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phase (Palla & Stahler 1993). This lack of knowledge is 
due to the fact that it is difficult to identify the earliest 
stages of the high-mass star formation process. Owing to 
their smaller numbers and shorter evolutionary timescales, 
high-mass star-forming regions are on the average more 
distant than low mass star-forming regions. Additionally, 
the clustered mode in which massive star formation gen- 
erally seems to proceed makes it very difficult to resolve 
and locate young high-mass protostars, which are deeply 
embedded and usually show little or no near-infrared emis- 
sion (see Menten & Reid 1995 for the case of the BN/KL 
region in Orion). 

Recently, much attention has been dedicated to stud- 
ies of massive star formation and significant progress is 
being made. A handful of sources likely to be young 
high-mass stellar objects have been identified. Among the 
best examples are G31. 41+0. 31mm (Cesaroni et al. 1994), 
G192.16 (Shepherd et al. 1998), IRAS 20126+4104 (Ce- 
saroni et al. 1997), IRAS 23385+6053 (Molinari et al. 
1998b) and G34.24+0.13mm (Hunter 1998). These ob- 
jects are characterized by high luminosities (> 10 4 L Q ), 
dense (> 10 6 cm -3 ) and warm (> 100 K) molecular gas 
and strong dust emission. Since they only show very weak 
or no free-free emission at cm wavelengths from ionized 
gas, they have not yet developed an ultracompact Hil 
(UCHii) region. Henceforth, we will call these objects in 
a pre-UCHu region phase high-mass proto-stellar objects, 
HMPOs for short. In this sense internally heated hot cores 
are a sub-class of HMPOs shortly before ultracompact Hn 
regions become detectable. 

The first cases of known HMPOs arc in the vicinity of 
UCHn regions, a selection effect because they were mostly 
discovered during molecular line and/or dust continuum 
studies of the UCHn regions. However, the nearness of an 
UCHn region (i.e. within, say, less than 0.1 pc or a few 
arcsec in angular units) is obviously a draw-back if one 
wants to study the earliest and unimpaired stages of mas- 
sive star formation: the physical and chemical conditions 
in the vicinity of the HMPOs are significantly modified by 
interactions with the UCHn regions and their precursors. 
Also, frequently it is obscrvationally difficult to discrimi- 
nate between material associated with the UCHn region 
and the nearby HMPO. The latter point can be partially 
resolved by high resolution and high dynamic range inter- 
ferometric observations, which in some of the cases cited 
above have led to the detection of HMPOs in the same 
fields as the UCHn regions. 

Systematic studies of a substantial sample of isolated 
HMPOs arc highly desirable and samples of HMPO can- 
didates have been investigated, e.g. by Molinari et al. 
(1996) and Hcnning et al. (2000). Here we describe a 
program of extensive observations of a new sample of can- 
didate HMPOs and present results indicating that some of 
the sources studied are objects in the earliest evolutionary 
stage known so far. Progress reports on this project were 
presented earlier by Menten et al. (1999), Sridharan et al. 
(1999) and Beuther et al. (2000). 

2. THE SAMPLE 

We have selected a sample of candidate HMPOs based 
on the criteria discussed by Ramesh & Sridharan (1997). 
They studied the reliability of the criteria defined by Wood 



& Churchwell (1989a) for identifying UCHn regions by 
their far-infrared properties and obtained a sample of rel- 
atively isolated HMPO candidates in the process. The 
objects are chosen from the IRAS Point Source Catalog 
and 

• are detected in the high density gas tracing CS 
J = 2 — > 1 survey of massive star forming regions 
(Bronfman et al. 1996). These sources satisfy the 
criteria used by Wood & Churchwell (1989) to select 
UCHn regions by their FIR colors. 

• are bright at FIR wavelengths (F eo > 90JyAFi O o > 
500Jy), 

• are not detected in the Galaxy-wide (single dish) 5 
GHz 1987 Green Bank (Gregory & Condon 1991) 
and Parkes-MIT-NRAO radio continuum surveys 
(Griffith et al. 1994; Wright et al. 1994) at flux 
densities above 25 mJy, 

• are north of —20° declination. 

The first two of these criteria are based on the expec- 
tation that dense gas is present and HMPOs and UCHn 
regions should resemble each other with regard to their 
dust and gas temperatures and luminosities. The third 
criterion ensures that we are dealing with isolated HMPO 
candidates in the sense that they do not have UCHn re- 
gions or more evolved Hil regions in their vicinity (within 
arcminutes, i.e. a few parsecs at the typical distances of 
a few kpc). Ramesh & Sridharan (1997) have shown that 
the surveys under consideration are sensitive enough to 
detect the radio continuum emission from typical UCHn 
regions throughout the Galaxy. The last criterium ensures 
the observability with the telescopes we are using. 

These criteria led to the identification of 69 candidate 
objects with kinematic distances (Brand & Blitz 1993) de- 
rived from the CS velocities. For sources inside the solar 
circle, there are two solutions for the kinematic distance. 
This ambiguity can in some cases be resolved by other 
means (see below). The sample is listed in Table 1 with 
systemic velocities, kinematic distances, and luminosities 
and temperatures derived from the HIRES database (see 
§3.1.1). Molinari et al. (1996) used a different set of far- 
infrared selection criteria to identify possible HMPO can- 
didates. Their and our sample have a total of 15 sources 
in common. 

Figure 1 (a and b) present plots of the Galactic distribu- 
tion and scale height for the whole sample. Sources with- 
out distance ambiguity are those in the outer Galaxy and 
those at the tangential points. Furthermore, we solved 
the distance ambiguity for a few more sources by other 
approaches: for sources with scale heights over 170 pc 
at the far distance [> 4 times the average Galactic scale 
height of UCHn regions, 37 pc (Bronfman et al. 2000)] the 
near distance was adopted (18151-1208, 18272-1217 and 
18517+0437); we chose this scaleheight, because 166 pc 
is the largest offset from the Galactic plane for a source 
with known distance (05553+1631). For 19217+1651 the 
far distance was adopted, because the observed VLA-cm 
fluxes are not explicable by an implied low mass star if the 
near distance is chosen. For a few objects distances of as- 
sociated sources were found in the literature (see Table 1). 
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We note that our selection criteria will certainly miss 
interesting classes of HMPOs. The first class we would 
miss are the HMPOs that are so young and cold that their 
SEDs peak longward of ~ 200 /im and are extremely weak 
at near infrared-wavelengths, making them undetectable 
by IRAS at 25 /im and thus excluding them from our 
sample. For example, a 10 4 M core of mean density of 
5 x 10 5 cm~ 3 and a temperature of 25 K, emitting as a 
blackbody at a distance of 5 kpc, produces a 25 /im flux 
density that is two orders of magnitude below the 0.5 Jy 
detection limit of IRAS at that wavelength. A prototype of 
such a cold, young object might be NGC 6334/I(NORTH) 
(Gezari 1982; Megeath & Tieftrunk 1999). The other 
missing group consists of sources with strongly associated, 
but not coincident, mid-infrared sources which increase the 
IRAS 12/zm flux and decrease the 25/xm/12/zm flux ratio, 
thus excluding them from our sample (e.g. 23385+6053, 
Molinari et al. 1998a). 

In this paper, the first of a series describing our HMPO 
survey, we give an overview of the overall characteristics 
of the sample, covering the mid-infrared (IRAS and MSX) 
and mm and cm continuum and spectral line observations. 
Spatial structures are discussed and we compare bolomet- 
ric luminosities, gas masses and cm luminosities of our 
sample with those of the ultracompact Hn regions. H 2 
and CH 3 OH maser observations as signposts of high-mass 
star formation are presented and we discuss spectral line 
observations of NH 3 , CO, CH 3 CN, thermal CH 3 OH, SiO 
and H 13 CO+. 

A more detailed analysis and discussion of the various 
datasets and follow-up observations on individual sources 
will be published in forthcoming papers. 

3. OBSERVATIONS AND RESULTS 
3.1. Archival data 
3.1.1. The IRAS-HIRES database 

The Infrared Astronomical Satellite (IRAS) performed 
an unbiased, sensitive all sky survey at 12, 25, 60 and 
100 /mi. Advanced processing of the survey data using the 
Maximum Correlation Method resulted in HIRES (High 
RESolution Processing) images, which have a higher spa- 
tial resolution than the original IRAS dataset (Aumann 
et al. 1990). We obtained images for the whole sample 
from the Infrared Processing and Analysis Center (IPAC) 
at 12, 25, 60, and 100 [im using 20 iterations of the algo- 
rithm. Spatial resolutions are between ss 40" (12 /im) and 
100" (100 /im). These are approximate values because the 
HIRES resolution varies around the sky. Detailed descrip- 
tions of HIRES can be found on the IPAC web-site 1 . 

In a majority of fields the 12 /im peak positions corre- 
spond well with the 100 /im peaks within the positional 
uncertainties (for an example see Fig. 2 bottom), but in 6 
of the sources the 12 /jm peaks are clearly offset by more 
than 1' from the 100 /im peaks (see Table 3 and example 
in Fig. 2 top). In some cases the 12 /im images are trac- 
ing more evolved IR sources in the vicinity of our targets 
(e.g. 05358+3543, Porras et al. 2000). In other cases the 
12 /urn peaks could be due to external heating of small 
grains (see discussion below). A number of sources are 

1 http://www.ipac.caltech.edu/ipac/iras/iras.html 

2 http://www.ipac.caltech.edu/ipac/msx/msx.html 



resolved at the higher resolution of the 12 /im images but 
appear unresolved at 100 /im (e.g. 05553+1631, see Fig. 2 
center). Thus, the fluxes in the different bands in the IRAS 
point source catalog do not necessarily correspond to the 
same objects. The real fraction of non-coinciding 100 and 
12 /xm sources might be even higher, because the spatial 
resolution, even of the HIRES images, is not very good. To 
improve on the Point Source Catalog data we obtained the 
fluxes around the main dust cores (Beuther et al. 2001) 
from the HIRES images and fitted two-component grey- 
bodies to these data (opacity corrected, modified black- 
bodies to fit also the 1.2 mm data, Beuther et al. 2001). 
A cold dust component fits the 60 /im and 100 /im data and 
a hot component fits the 12 /im and 25 /im data. These 
fits reproduce the HIRES data well (see Figure 3 for ex- 
amples) for temperatures of the cold dust component T c( j 
around 50 K and hot component Thd between 150 K and 
200 K (Table 1). While the hot component might repre- 
sent a more evolved and hot inner core, it could also be 
due to small non-thermally excited grains exposed to ex- 
ternal UV heating as found for many other IRAS sources 
(Mathis 1990). The MSX data (see §4.1) support this in- 
terpretation for a fraction of sources for which the MSX- 
and mm positions do not coincide. Integrating the two- 
component fits gives estimates of the total luminosities of 
our objects (Tabic 1). We find total luminosities in the 
range 10 35 — 10 6 L Q , therefore we conclude that the cm- 
bedded sources are high-mass objects with ZAMS masses 
> 8M Q (spectral type earlier than B2). The ratio of the 
flux in the hot component to that in the cold component 
varies between 0.05 and 0.4 with an average value around 
0.16. This shows that the main contribution to the total 
luminosity comes from the cold components. 

3.1.2. The Midcourse Space Experiment (MSX) Point 
Source Catalog 

MSX, a Ballistic Missile Defense Organization satellite, 
surveyed the whole Galactic plane (b±6°) in 6 bands cen- 
tered at 4.29, 4.35, 8.28, 12.13, 14.65, and 21.34 /im. 
The sensitivity in the 4 /im bands is between 10 and 
30 Jy/beam, in the other bands the sensitivity is better, 
ranging between 0.1 and a few Jy/bcam. With a spatial 
resolution of 18.3", positions in the Point Source Catalog 
are accurate to an rms of 4" — 5", and the calibration er- 
rors are well within 13%. For a more detailed description 
see, e.g., Egan et al. (1998) and the MSX web -site 2 . 

The MSX Point Source Catalog lists sources in 62 out 
of our 69 fields (Table 1). In a few cases multiple sources 
are found and a total of 87 sources are found toward our 
target regions. Except for 18540+0220 and 23151+5912 
no sources were detected in the 4 /im bands, while in the 
other bands almost all of the sources are detected. Com- 
paring the fluxes in the 12 /im and 21 /im bands (nearest 
to the 12/im & 25/im IRAS bands) the flux ratio is be- 
tween 0.1 and 0.5 in 63 MSX sources, between 0.6 and 0.8 
in 4 sources and undefined in 20 because of non-detection 
in either of the two bands. The flux ratios between 0.1 
and 0.5 are similar to the IRAS flux ratios. For the main 
cores the flux ratios do not vary between sources with or 
without cm detections, therefore the MSX fluxes alone do 
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not allow determinations of the evolutionary phase of a 
source. The four sources with ratios in excess of 0.5 are 
found at the edge or offset from the main core as seen in 
the mm continuum and are most likely hotter and more 
evolved sub-sources. 

3.2. New observations 

3.2.1. 1.2 mm dust continuum emission 

Dust emission produced by massive cores has been found 
and mapped in all the fields and a detailed analysis of this 
extensive database is presented in an accompanying paper 
(Bcuther et al. 2001). For the present discussion we only 
use the core masses derived from the 1.2 mm data. 

3.2.2. 3.6 cm continuum emission 

Because small Hn regions are in general difficult to de- 
tect in most existing centimeter-wavelength single dish 
surveys of the Galactic plane, we used the Very Large Ar- 
ray (VLA) 3 to search for weak free-free emission at 3.6 cm 
using the B array with a synthesized beam of 0.7" and a lcr 
rms sensitivity of ~ 0.1 mJy. The observations were car- 
ried out on July 2, 1998. Each source was observed twice 
at widely different hour angles, for a total of 10 minutes in 
order to improve uv-coverage. Suitable phase calibrators 
were observed at intervals of approximately 30 minutes. 
The data were reduced using NRAO's AIPS software us- 
ing standard procedures. 

Radio emission from objects in the pre-UCHu phase 
might be caused by recently ignited massive protostars 
and/or jet phenomena (e.g. Hofner et al. 1999). Back- 
ground contamination by extragalactic radio sources at 
a 1 mJy level is estimated to be extremely low at this 
frequencies (0.007 per arcmin 2 , Fomalont ct al. 1991). 
Other processes such as synchrotron emission (Reid et al. 
1995) or stellar winds could be producing the observed 
emission which require spectral index information for fur- 
ther elucidation. About 10% of the sources originally ob- 
served by Wood & Churchwell (1989b) are optically thick 
at 2 cm, and a similar or even higher percentage - since in 
the early stages we expect to find hypercompact sources 
(< 2000 AU, Tieftrunk et al. 1997) - might be applica- 
ble for our detections. Sources without any cm emission 
could be even younger objects, where either no star has 
ignited or the ionized region is so young and small that 
the free-free emission is still confined and not detectable 
down to the mJy level (Churchwell 2000). The mm ob- 
servations (Beuther et al. 2001) prove that our sample in- 
deed consists of massive cores which should produce high- 
mass stars, assuming normal IMFs. Out of the 66 fields 
observed at 3.6 cm, 24 show no emission at the 1 mJy 
level, and 19 have weak emission with flux densities below 
10 mJy, while 21 sources show emission between 10 and 
a few 100 mJy. For the remaining 2 fields the dynamic 
range of our snapshot data is not good enough to draw a 
firm conclusion. Observed flux densities are listed in Tabic 
3. Typical ultracompact Hn region fluxes at distances of a 
few kpc are between a few mJy and more than a Jy (Wood 
& Churchwell 1989b), which is far higher than the flux 
densities observed in a substantial fraction of our sample. 
But there are also a number of well known sources which 

3 The VLA is operated by the National Radio Astronomy Observatory 
with the National Science Foundation (NSF). 



exhibit cm continuum emission at a level not detectable 
with our sensitivity. The W3(OH) Turner- Welch object 
(also known as W3(H 2 0)) at a distance of 2.2 kpc has a 
cm flux of 1.1 mJy (Reid et al. 1995) due to synchrotron 
emission, which is similar to our detection limit, but other 
sources like source I in Orion KL (Menten & Reid 1995) or 
Cep A East (Garay et al. 1996) would show fluxes below 
0.3 mJy if moved to a distance of 2 kpc, which is below our 
detection limit. Therefore it is likely that deeper images of 
our sample would reveal more cm sources. The fact that 
the sample contains sources at different stages of evolu- 
tion makes it an ideal laboratory to explore the evolution 
of massive stars with a consistent database. 

3.2.3. Spectral line observations 

In order to understand the chemistry, kinematics and 
energetics of our target regions, we conducted a number 
of different spectral line observations from mm to cm wave- 
lengths. 



Millimeter observations: outflow and dense core tracers 

The IRAM 30-m telescope was used (mainly in spring 
1998, and for several periods of a few hours duration in 
spring 1999 and 2000) to obtain single point spectra for 
36 sources in a number of tracers sensitive to the physi- 
cal and chemical conditions expected near HMPOs: high 
temperatures and densities, elevated abundances of molec- 
ular species resulting from evaporation of icy grain man- 
tles, and outflow activity. In most cases a number of lines 
could be observed simultaneously (see Table 2). 

The observations were conducted in the wobbler switch- 
ing mode, where the wobbling secondary reflector switched 
in azimuth between the ON and an OFF position with a 
switching amplitude of 180" and a frequency of 0.5 Hz. For 
the CO data the OFF position is generally not free of emis- 
sion, which corrupts the spectra at velocities around the 
systemic velocity but in general does not affect spatially 
more confined broad line wings in which wc are primarily 
interested. As backends, we used the facility autocorrela- 
tor and two 512 MHz filtcrbanks. 

To complete the CO line wing/outflow single positions 
search we observed the rest of the sample in June 2000 at 
the CSO in the 12 CO(2 1) line (see Table 2). These 
data were obtained in the position switching mode, where 
the OFF positions were determined using Galactic CO sur- 
veys (Sanders et al. 1986; Dame et al. 1987) and verified 
to be free of emission. 

We found CO line wings tracing bipolar outflows in 58 
out of 69 sources. Six sources showed no line wings, and for 
5 sources the spectra are too confused by Galactic plane 
emission to allow a classification. Assuming that line wings 
are only detectable for outflow sources with an inclination 
angle of at least 10° to the plane of the sky, statistically 
11% (in our case « 8 sources) should show no wing emis- 
sion although an outflow could be present. That would be 
consistent with outflows being present in all of our sources 
and strongly suggests that outflow activity is an ubiquitous 
phenomenon not only in low mass star-forming regions but 

AO) with Associated Universities, Inc., under cooperative agreement 
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also in regions of massive star formation. Previous inves- 
tigations that obtained similar results are e.g. Shepherd 
et al. (1996); Henning et al. (2000); Zhang et al. (2001) 
(see also the review by Churchwell 2000). 

SiO emission as another outflow tracer (e.g. Schilke et 
al. 1997) was detected in the 30-m observations toward 23 
out of 36 sources with wing emission seen in most cases. 

Thermal emission from CH 3 OH and/or CH 3 CN was de- 
tected in 19 out of the 30 sources implying that hot core 
activity is prevalent. H 13 CO + (l — ► 0) emission was found 
in all 36 observed sources, most likely tracing structures 
similar to the dust cores. We mostly used IRAS coordi- 
nates as source locations because at the time of the ob- 
servations improved positions from 1.2 mm observations 
were not yet available. In 8 fields the main mm peak is 
offset from the IRAS position by more than the 1.2 mm 
FWHM beamwidth of 11", which would have lowered the 
detection rates. 

H 2 and CH^OH masers 

During 1998 February and November, we searched for 
maser emission in the 22 GHz H2O (mean rms « 0.4 Jy) 
and 6.7 GHz CH3OH lines (mean rms ~ 0.2 Jy) toward 
all of our target sources using the Effelsberg 100 m tele- 
scope (see Table 2). The calibration error is estimated to 
be around 20%. 

Both maser searches had high detection rates: 29 H2O 
(10 of them new detections) and 26 CH3OH masers (5 new 
detections) were found; toward 19 sources both H2O and 
CH3OH masers were detected (see Fig. 5 for examples). 
The CH3OH detection rate of our sample differs from that 
of UCH11 regions: We find CH 3 OH maser toward 38% of 
the sources in our sample while Walsh et al. (1998) quote 
a detection rate of 20% toward UCH11 regions. In contrast, 
the detection rate of the H 2 masers toward the sources 
of our sample is only 42%, while Churchwell et al. (1990) 
detect H2O masers toward 67% of 84 UCH11 regions they 
surveyed with comparable flux limits. 

In Figure 6, we present the source and maser distri- 
bution of our sample in a Wood & Churchwell (1989)- 
style color-color diagram. While H2O and CH3OH maser 
sources are found all over the diagram, combined clus- 
tering of both maser types is much more common in the 
top-right corner, corresponding to the coldest and hence 
presumably the least evolved sources, which makes sources 
with colors in that area the most promising HMPO can- 
didates. 

While methanol maser emission is always confined to 
less than 15 km s _1 in total around the systemic velocity 
the water masers nearly always emit over a much wider 
range (up to 70 km s _1 ). Just comparing the same ve- 
locity ranges both types mostly show emission features 
at different velocities, which strongly confirms earlier as- 
sumptions that these maser types are emitted in different 
regions of the star-forming core (Menten 1996). While the 
methanol masers are probably confined to the innermost 
part - the molecular envelope of a central star or possibly 
a disk - water maser emission could originate there as well, 
but the higher velocities point to different emission spots, 
most likely in outflows (e.g. Marti et al. 1999). 

Follow-up interferometric observations with the VLA 
(H 2 0) and the Australian Telescope Compact Array 



(ATCA, CH3OH) have been conducted and a detailed 
analysis and comparison of the data with the other 
datasets will be published soon. 

NH3 observations 

To obtain estimates of the kinetic temperature, the 100- 
m telescope was used to observe the NH 3 (J, K) = (1,1), 
(2,2) inversion lines in the entire sample (sec Table 2). 
The (1,1) line was detected toward 59 sources, and for 
40 sources the (2,2) line was found to be strong enough 
to allow a meaningful estimate of the gas temperature. 
The rotation temperatures listed in Table 1 were derived 
following Ungerechts et al. (1986); in the relevant tem- 
perature range around 20 K kinetic temperatures are only 
marginally higher (Danby et al. 1988). The average ro- 
tation temperature is 19 K with the notable exception of 
18089-1732 for which we derive 38 K. Molinari et al. (1996) 
found a similar mean value of T^ n ~ 22 K in their sample 
of young massive objects, which is identical to the aver- 
age rotation temperature Churchwell et al. (1990) found 
toward UCH11 regions. Comparing the NH 3 derived tem- 
peratures with the dust temperatures (from HIRES, Table 
1) it is evident that the derived dust temperatures are on 
the average ~ 20 K higher than the rotation temperatures. 
This difference reflects the fact that IRAS was not sensitive 
to temperatures lower than 30 K and therefore detected 
emission from warm dust inside the core. In contrast to 
that the (1,1) and (2,2) ammonia lines trace the cooler 
and more extended envelope and therefore result in lower 
temperatures. It is possible and likely that NH3(4,4) ob- 
servations will reveal higher temperatures due to a warm 
inner core. 

The average intrinsic NH3 line widths we find in our 
sample is 2.0 km s _1 , while the average value found toward 
a sample of molecular cores associated with UCH11 regions 
is 3.1 km s _1 (Churchwell et al. 1990). This indicates sig- 
nificantly less turbulence in the presumably younger HM- 
POs possibly because in this young sources there has been 
less input of mechanical energy yet, e.g. by outflows. 

4. DISCUSSION 

4.1. Source Multiplicity 

Positional comparison between MSX-/mm sources 
(Beuther et al. 2001) and MSX-/cm sources (Table 3) 
reveal median separations of approximately 7" (between 
1" and 30") in both cases. The positional accuracy of 
MSX is 4" - 5" (Egan et al. 1998), comparable to that of 
the mm data, while the the cm data are more accurate, at 
better than 1". While the offsets mentioned might be due 
to positional uncertainties in many cases, there remain sig- 
nificant numbers of fields where the offsets are larger than 
any observational error and therefore the different tracers 
depict different sources. MSX sources located at the edges 
of the cores observed at mm-wavelengths might either be 
reflection nebulae, i.e. reprocessed radiation from a deeply 
embedded source, or they could also be sources belonging 
to the same evolving cluster as the mm source, but in a 
more advanced stage of evolution. Even younger stages of 
evolution might even be undetectable at 12/xm (see also 
§2). 
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In 29 out of 42 fields with cm detections, MSX-, mm- 
and cm-sources are found. In most of these cases it is dif- 
ficult to draw firm conclusions as to whether the emission 
in the three wavelength regimes arises from the same or 
from different sources, because the spatial resolution is in- 
sufficient. For the regions showing extended cm emission 
(of size larger than a few arcsecs), in 5 out of 7 sources the 
MSX-source is found to be associated with the cm- rather 
than with the mm-source. While we are dealing with low- 
number statistics, different stages of evolution in the same 
cluster are the most probable cause. The more evolved and 
therefore hotter sources exciting an UCHn region are more 
likely to emit at 12 ^.m and thus be detectable by MSX 
than the colder and more deeply embedded mm sources. 

Generally, sources showing cm emission are found to be 
associated with the mm cores (Beuther et al. 2001), but 
the cm peak is mostly offset from the mm peak (Table 3). 
The median offset between mm and cm sources is 11" (be- 
tween 2" and 40"). Given the positional accuracies quoted 
above, these offsets are real and not due to pointing prob- 
lems. Assuming the near distance for sources with distance 
ambiguities the median projected separation between mm 
and cm peaks is « 0.1 pc. Thus, while the sources most 
likely belong to the same (proto)cluster, they are clearly 
different sub-sources. The cm sources are probably slightly 
older and more evolved, while the mm sources might rep- 
resent the youngest objects or object groups in the fields. 
For the few sources which seem to be coincident possible 
projection effects have also to be considered. 

4.2. Mass-luminosity relations 

To understand mass and luminosity evolution in the 
massive star formation process, we compare the masses de- 
rived from the mm continuum data (Beuther et al. 2001) 
and luminosities from the HIRES-data of our sample with 
other classes of objects. Since we selected our sample to 
contain objects in a pre-UCHu region phase, natural ref- 
erence objects are UCHn regions. For this purpose we use 
UCHn region data from Hunter (1997) and Hunter et al. 
(2000) and estimate masses from 350/im images using the 
same assumptions as employed by us (Beuther et al. 2001) 
for our 1.2 mm data. For a few sources Hunter et al. also 
have 1.3 mm data, and the masses separately derived us- 
ing only the 1.3mm fluxes and the 350^m fluxes agree well 
with each other. 

In Figure 4(a) we plot far-infrared luminosities versus 
core masses for HMPOs and UCHn regions. The two 
groups appear separated, with the UCHn regions hav- 
ing higher bolometric luminosities than the HMPOs for 
the same core masses. We interpret this as an evolution- 
ary effect caused by the embedded cluster, which becomes 
more luminous at FIR wavelengths while they evolve and 
destroy the cores from which they were born. Another 
possible explanation for this trend would be that different 
massive star clusters differ in their initial mass function 
(IMF) and that sources with UCHn regions form more 
massive stars than HMPOs. We think this is unlikely, 
because Massey (1998) finds for more than 20 OB associ- 
ations essentially a Salpeter IMF (N cx to -135 ) with little 
variation. Therefore, we argue that we are indeed observ- 
ing an evolutionary effect and our sample is on the average 
younger than UCHn regions. 



Monte Carlo simulations to derive expected cluster lumi- 
nosities over a wide range of cluster masses were performed 
using the methods outlined in Walsh et al. (2001). As Ini- 
tial Mass Functions we used a Salpeter IMF (N cx m -135 ) 
for sources above 1 M Q (Massey 1998) and a Scalo IMF 
(N cx to~ 0,83 ) for sources below 1 M . The upper end of 
the IMF is still under discussion and Casassus et al. (2000) 
propose a steepening there, but we work with the distri- 
bution found by Massey (1998), which is based on obser- 
vations of more than 20 OB associations. A steeper IMF 
would flatten the top end of the derived cluster luminosity 
distribution slightly, but qualitatively the results stay simi- 
lar. As mass-to-luminosity relation L cx vnP we use (3 = 2.8 
for m < 1, /3 = 4 for 1 < m < 30 and (3 — 2 for m > 30 
(Massey 1998; Schatzman & Praderie 1993). To scale 
the cluster masses to the initial core masses we assumed 
a star formation efficiency of 30% (Lada 1993), and the 
results are shown in greyscale in Figure 4(a). The shaded 
area indicates where 90% of the star clusters are expected 
to be found. The large spread in luminosity with given 
mass can be explained by the strong power law depen- 
dence between mass and luminosity. Many low mass stars 
do not raise the luminosity of the the cluster much, but the 
formation of just 1 massive star increases the cluster lumi- 
nosity over orders of magnitude. The upper cutoff at lumi- 
nosities higher than 10 6 L Q occurs because of an artificial 
upper mass truncation in the simulated IMF at 100 M . 
The main conclusion we derive from these simulations is 
that many HMPOs arc undcrluminous compared to the 
expected cluster luminosities. The UCHn region luminosi- 
ties agree better, but the most massive UCHn regions are 
found to be undcrluminous again. Since we expect the 
same formation processes in both types of sources with 
the same underlying IMFs (as stated by Massey 1998), 
this strengthens our conclusion that the HMPOs are in a 
younger stage of evolution and have not yet reached their 
final cluster luminosity. The Figure suggests further that 
objects evolve from HMPOs to UCHn regions to clusters. 

A further support of the youth of our sample is outlined 
in Figure 4(b), which presents the distance independent 
M/L ratio in a histogram. As already expected from Fig- 
ure 4(a) the ratio M/L is lower on the average for UCHns 
(peak at 0.01) than for HMPOs (peak at 0.05). This 
again can be explained with the HMPOs being younger 
and therefore less luminous than the more evolved UCHn 
regions. A higher M/L ratio could also indicate sample of 
less massive stellar objects on an average, but this is not 
the case in our sample because the core mass ranges of the 
HMPOs and the UCHn regions are similar (see Fig. 4(a)). 

4.3. Bolometric- versus cm luminosities 

We calculate from the measured cm fluxes the expected 
stellar luminosities of the central sources [as outlined in 
Kurtz et al. (1994, cqs. 1 & 3) assuming the emission to 
be optically thin at 3.6 cm] and compare them to the IR 
luminosities obtained from the IRAS data base. For opti- 
cally thick emission, which we expect in a fair number of 
our sources, the derived luminosities are lower limits. Fig- 
ure 7 presents the derived cm luminosities plotted against 
the corresponding IR luminosities. It is obvious that the 
IR luminosities are on the average far higher than the cm 
luminosities; similar results were found by Testi (2000) 
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and Walsh et al. (2001) on different samples. This most 
likely reflects the fact that the cm data trace just one sin- 
gle already ignited star while the IR data - due to the 
large IRAS beam - are tracing a whole cluster with many 
sources, some younger with no H-burning at their center 
yet, and others too weak to ionize a significant UCHn re- 
gion, which by now is the standard scenario. This would 
be consistent with massive stars forming only in clusters 
(Stahler et al. 2000). A possible caveat is that some of 
the young ultracompact or even hypercompact Hn regions 
could be confined strongly and be optically thick which 
leads to an underestimate of the cm luminosities. Data 
at other wavelengths are necessary to define the spectral 
indices and by that the ratio of optically thick sources. Ad- 
ditionally, accretion contributions to the luminosities have 
to be taken into account (see §4.4). All these phenomena 
explain the observations qualitatively well, but separating 
the different contributions quantitatively is not possible 
with the present database. 

4.4. Accretion luminosity 

Luminosity due to accretion plays a significant role in 
the formation of massive stars. We assume a simplified 
relation for the accretion luminosity (Wolfire & Cassinelli 
1987): 

_ GM*M 

-^acc D 

with G the gravitational constant, M* the mass of the 
(proto)star, M the accretion rate and i?* the stellar radius. 
Main sources of uncertainty are the accretion rate and the 
radius of the accreting star. Assuming an accretion rate 
of M = 10~ 3 M Q /yr (based on outflow observations of 
this sample presented in a forthcoming paper by Beuthcr 
et al., in prep.) the accretion luminosity for stars between 
10 and 40 M Q is around 6 x 10 4 L . While the accretion 
rate we choose for this estimate is an upper limit in cur- 
rent theories of high-mass star formation (e.g. Bernasconi 
& Maeder 1996; Norberg & Maeder 2000; Osorio et al. 
1999), the radius of the evolving protostar is expected to 
be larger than the final main sequence radius we employed 
(Bernasconi & Maeder 1996). Thus, these accretion lu- 
minosities are upper limits to true accretion luminosities. 
With this assumed accretion rates it takes 3 x 10 4 years 
to build a 30 M Q star which has a luminosity comparable 
to the accretion luminosity. Using the increasing accre- 
tion rates by Norberg & Maeder (2000) equality of stellar 
and accretion luminosities is reached even later, at times 
around 10 6 years. Osorio et al. (1999) find for a num- 
ber of well known hot cores that accretion still produces 
around 90% of the bolometric luminosity after a few times 
10 4 years. In the models of intermediate-mass stars by 
Palla & Stahler (1993) the radiative luminosity matches 
the accretion luminosity after around 10 5 years. While 
the time scales differ between the different models, in all 
the cases the accretion luminosity plays a dominant role 
for a significant amount of time during the massive cluster 
forming process. Now, comparing the accretion luminosi- 
ties with the IR derived luminosities of our sample (Table 
1) it seems very likely that in many objects a significant 
part of the detected flux is due to accretion. However, at 
this stage of the investigation we cannot determine what 



fractions of the total luminosity are due to accretion and 
nuclear burning. High accretion rates (Norberg & Maeder 
2000) would also cause quenching of the UCHn cm radi- 
ation in the initial stages of evolution (Walmsley 1995; 
Osorio et al. 1999), consistent with our low average cm 
fluxes. 

5. CONCLUSIONS 

Our selection criteria - basically CS detections of bright 
IRAS point sources with colors similar to ultracompact 
Hn regions and the absence of significant free-free emis- 
sion at 5 GHz - were very successful in finding sources 
at the earliest known stages of high-mass star formation. 
An analysis of the JiMS-HIRES and MSX-data gives an 
overview of spatial structures on large scales and the spec- 
tral energy distributions. We find bolometric luminosities 
of 10 4 — 10 6 L Q and dust temperatures around 40 K. 

• VLA cm observations reveal that a large fraction of 
the sample does not emit free-free emission down 
to 1 mJy, which suggests that the sources are at 
very early evolutionary stages. The majority of 
the detected cm peaks is spatially offset from the 
mm peaks which implies different sources in dif- 
ferent evolutionary stages being responsible for the 
different emission features. 

• Luminosities based on the cm emission are on the 
average far lower than bolometric luminosities de- 
rived in the FIR. The most likely explanation for 
this phenomenon is that the cm emission traces just 
one massive star whereas the FIR-emission reflects 
the whole cluster luminosity. It has to be taken into 
account that the cm emission might still be opti- 
cally thick in a number of sources, which results in 
an underestimation of the intrinsic luminosity. 

• HMPOs are found to have a higher mass to luminos- 
ity ratio than UCHn regions. The most likely expla- 
nation is that the HMPOs are significantly younger 
with more molecular gas left and less (stellar) ra- 
diation reprocessed in the IR. Evolving further to 
the UCHn region stage the sources become more lu- 
minous and start destroying their surrounding gas 
cocoons. 

• Comparing estimated accretion luminosities with 
the derived bolometric luminosites indicates that in 
many sources a large fraction of the total luminosity 
could be due to accretion processes. 

• Outflow phenomena are ubiquitous and methanol 
and water maser detection rates are high. Other 
molecular species such as SiO (outflow tracer), 
H 13 CO + (dense gas tracer) and hot core tracers like 
CH 3 OH and CH3CN were found in a large fraction 
of the observed sources. 

• Rotation temperatures from NH3 (1,1) and (2,2) 
lines in 40 out of 69 sources are on the average 20 K, 
about 20 K lower than the HIRES-derived dust tem- 
peratures. This phenomenon has been found fre- 
quently and comes about because IRAS was not 
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sensitive to temperatures below 30 K and observed 
warmer dust further inside the core. On the con- 
trary, NH 3 (1,1) and (2,2) is sensitive to the cooler 
and more extended envelope, which results in lower 
temperatures. 

We have carried out further single dish and intcrferomet- 
ric studies of promising objects from this sample. These 
results will be published in forthcoming papers (Bcuthcr 
et al. Wyrowski et al. in preparation). Combined with ef- 
forts by other groups a reasonable sized sample of HMPO 
candidates now exists, and in the near future a good ob- 
servational characterization of HMPOs can be achieved to 
confront theories. 

We are indebted to Andrew Walsh for the Monte Carlo 
simulations of the cluster luminosities. We thank Frank 



Bertoldi, Malcolm Walmslcy and Floris v.d. Tak for many 
useful comments on earlier versions of this paper. Addi- 
tionally, we thank Eberhard Hansis who worked on the 
HIRES data and resolved some distance ambiguities. We 
also like to thank Frederique Motte, Bernd Weferling and 
again Frank Bertoldi, who conducted the bolometer ob- 
servations at Pico Veleta. Last but not least, we thank 
the referee for additional improving comments. The ob- 
servational work was conducted at the Effelsberg 100m, 
the VLA, the IRAM 30m and the CSO (suppported by 
the NSF grant AST 96-15025). It has made extensive use 
of the IPAC and SIMBAD databases. H. Beuther gets 
support by the Deutsche Forschungsgemeinschaft, DFG 
project number SPP 471. F. Wyrowski is supported by 
the National Science Foundation under Grant No. AST- 
9981289. 
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-08 


23 


35 


84.4 


10.4 


5.2 


40 
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4 


.8 


4.2 




18 


3 


18337 — 0743 


18 


36 


29.0 


-07 


40 


33 


57.9 


11.5 


4.0 


37 


139 


5 





4.0 




17 


1 


18345-0641 


18 


37 


16.8 


-06 


38 


32 


95.9 


9.5 5 




36 


161 


4 


6 






16 


1 


18348 — 0616 


18 


37 


29.0 


-06 


14 


15 


109.5 


9.0 


6.3 


39 


159 


5 


.1 


4.8 




19 


3 


18372 — 0541 


18 


39 


56.0 


-05 


38 


49 


23.6 


13.4 


1.8 


33 


JiU 




3 


3.5 


110 


16 


1 


18385 — 0512 


18 


41 


12.0 


-05 


09 


06 




13. 1 


2.0 




onl 


5 


.3 


3.7 








1 8-1.20 - 020 '1 




45 


12.8 


-02 


01 


12 




1 q Q 


« \ 


q5i 




5 





2-8 








1813 1 — 03 12 


i » 


45 


46.9 


-03 


09 


24 


OK o 




6.7 




1 oo 


4 


.5 


4.4 




15 


2 


18137 — 0216 


18 


46 


22.7 


-02 


13 


24 


lino 
110.8 






31 


221 


4 


.4 










18 1 10 — 0148 


18 


46 


36.3 


-01 


45 


23 


97.6 


6 

8.3 




97 


121 


4 


7 




75 


23 


1 


18415 — 0222 




47 


10.8 


-02 


19 


06 












5 


.2 










18447 — 0229 


18 


47 


23.7 


-02 


25 


55 


102 6 


8 2 


6 6 


36 


172 


4 


6 


4 4 




15 




18454 — 0136 


18 


48 


03.7 


— 01 


33 


23 


38.9 


11.9 


2.7 


31 


142 


4 


.8 


3.5 




22 


1 


18154 — 0158 




48 


01.3 


— 01 


54 


49 




5 6 7 








4 


.3 










18460 — 0307 


18 


48 


39.2 


-03 


03 


53 


837 




5 2 


43 


176 


4 


.6 


4 1 




19 


^ 


18470-0044 


18 


49 


36.7 


-00 


41 


05 


96.5 


8.2 7 




56 


158 


4 


.9 






20 


2 


18472-0022 


18 


49 


50.7 


-00 


19 


09 


49.0 


11.1 


3.2 


49 


135 


4 


9 


3.8 






1 


18488+0000 


18 


51 


24.8 


+ 00 


04 


19 


82.7 


8.9 


5.4 


42 


169 


4 


9 


4.5 




20 


1 


18517+0437 


18 


54 


13.8 


+ 04 


41 


32 


43.9 


2.9 




38 


176 


4 


.1 








1 


18521+0134 


18 


54 


40.8 


+ 01 


38 


02 


76.0 


9.0 


5.0 


37 


151 


4 


6 


4.1 


90 




1 


18530+0215 


18 


55 


34.2 


+ 02 


19 


08 


77.7 


8.7 


5.1 


51 


127 


5 


.4 


4.9 




16 


1 


18540+0220 


18 


56 


35.6 


+ 02 


2d 


54 


49.6 


10.6 


3.3 


38 


159 


4 


9 


3.9 






2 


18553+0414 


18 


57 


52.9 


+ 04 


18 


06 


10.0 


12.9 


0.6 


42 


152 


5 


.1 


2.4 








18566+0408 


18 


59 


09.9 


+ 04 


12 


14 


85.2 


6.7 




51 


161 


4 


.8 






15 


1 


19012+0536 


19 


03 


45.1 


+ 05 


40 


40 


65.8 


8.6 


4.6 


43 


139 


4 


.7 


4.2 




21 


1 


19035+0641 


19 


06 


01.1 


+ 06 


46 


35 


32.4 


2.2 8 




51 


100 


3 


9 






21 


1 


19074+0752 


19 


09 


53.3 


+ 07 


57 


22 


54.8 


8.9 


3.7 


45 


144 


4 


.8 


4.0 




16 


1 


19175 + 1357 


19 


19 


49.1 


+ 14 


02 


46 


14.6 


10. 6 9 




36 


165 


4 


8 








3 


19217+1651 


19 


23 


58.8 


+ 16 


57 


37 


3.5 


10.5 




38 


124 


4 


9 






25 




19220 + 1432 


19 


24 


19.7 


+ 14 


38 


03 


68.8 


5.5 




43 


151 


4 


1 






23 


1 


19266 + 1745 


19 


28 


54.0 


+ 17 


51 


56 


5.0 


10.0 


0.3 


32 


176 


4 


.7 


1.7 






2 


19282 + 1814 


19 


30 


28.1 


+ 18 


20 


53 


23.6 


8.2 


1.9 


36 


231 


4 


9 


3.6 






1 


19403 + 2258 


19 


42 


27.2 


+ 23 


05 


12 


26.7 


6.3 


2.4 


57 


150 


4 


.7 


3.8 






1 


19410 + 2336 


19 


43 


11.4 


+ 23 


44 


06 


22.4 


6.4 


2.1 


46 


147 







4.0 




18 


1 


19411 + 2306 


19 


43 


18.1 


+ 23 


13 


59 


29.0 


5.8 


2.9 


41 


157 


4 


.3 


3.7 




14 


1 


19413+2332 


19 


43 


28.9 


+ 23 


40 


04 


20.8 


6.8 


1.8 


41 


133 


4 


.4 


3.3 




18 


1 


19471 + 2641 


19 


49 


09.9 


+ 26 


48 


52 


21.0 


2.4 10 




38 


194 


3 


6 




130 




2 


20051 + 3435 


20 


07 


03.8 


+ 34 


44 


35 


11.6 


3.7 


1.6 


49 


143 


4 





3.3 






1 


20081 + 2720 


20 


10 


11.5 


+ 27 


29 


06 


5.7 


0.7 2 




34 


248 


2 


.5 








1 


20126+4104 


20 


14 


26.0 


+ 41 


13 


32 


-3.8 


1.7 2 




62 


112 


3 


9 






23 


1 


20205 + 3948 


20 


22 


21.9 


+ 39 


58 


05 


-1.7 


4.5 




48 


175 


4 


.5 








1 


20216+4107 


20 


23 


23.8 


+ 41 


17 


40 


-2.0 


1.7 2 




46 


159 


3 


.3 






21 


1 


20293 + 3952 


20 


31 


10.7 


+ 40 


03 


10 


6.3 


2.0 


1.3 


56 


157 


3 


.8 


3.4 




15 


1 


20319 + 3958 


20 


33 


49.3 


+ 40 


08 


45 


8.8 


1.6 




73 


162 


3 


.8 








1 


20332+4124 


20 


35 


00.5 


+ 41 


34 


48 


-2.0 


3.9 




56 


148 


4 


.4 






17 




20343+4129 


20 


36 


07.1 


+ 41 


40 


01 


11.5 


1.4 




44 


150 


3 


5 






18 




22134 + 5834 


22 


15 


09.1 


+ 58 


49 


09 


-18.3 


2.6 




61 


135 


4 


.1 






18 




22551+6221 


22 


57 


05.2 


+ 62 


37 


44 


-13.4 


0.7 3 




45 


186 


3 


.2 










22570 + 5912 


22 


59 


06.5 


+ 59 


28 


28 


-46.7 


5.1 




54 


145 


4 


7 








1 


23033+5951 


23 


05 


25.7 


+ 60 


08 


08 


-53.1 


3.5 4 




52 


149 


4 


.0 






20 


1 


23139 + 5939 


23 


16 


09.3 


+ 59 


55 


23 


-44.7 


4.8 




41 


140 


4 


.4 








1 


23151 + 5912 


23 


17 


21.0 


+ 59 


28 


49 


-54.4 


5.7 




68 


175 


5 


.0 








1 


23545+6508 


23 


57 


05.2 


+ 65 


25 


11 


-18.4 


0.8 2 




50 


181 


3 


.0 






18 


1 



Note. — Only far and no near distance means that the distance ambiguity is solved. Rotation temperatures T ro ^(NH3) are correct to approximately 3 K. The 12/100 column 
refers to offsets between the corresponding IRAS sources and the last columns shows the number of MSX point source catalog detections in fields of radius 2 around the IRAS 
position. 1 Snell et al. 1990; 2 Wilking ct al. 1989; 3 Blaauw 1964; 4 .Tijina ct al. 1999; 5 Garay ct al. 1993; 6 Walsh et al. 1997; 7 Kuchar & Bania 1994; 8 Ostcrloh ct al. 
1996; 9 near star HD231160; 10 near star CGO 544 
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Table 2 
Observing parameters 





trcq. 


HPBW 


'isys 


An 


Obs. 


tracer 


12 CO 2 1 


[GHz] 


["] 


[K] 


[km s" 1 ] 






230.5 


11 


250 


0.1 


PV 


outflow 


SiO 2 -. 1 


86.9 


29 


85 


3 


PV 


outflow 


H 13 CO + 1 -> 


86.9 


29 


85 


3 


PV 


dense core 


13 CO 1 -> 


110.2 


22 


120 


0.8 


PV 


dense core 


CH 3 OH J fc = 5 fc -> 4 fc 


~ 241.8 


11 


250 


0.1 


PV 


dense core 


CH 3 CN ,7 = 6^5 


~ 110.4 


22 


120 


0.8 


PV 


dense core 


12 CO 2 1 


230.5 


27 


500 


0.06 


CSO 


outflow 


NH 3 (1,1), (2,2) 


~ 22.6 


40 


50 


0.25 


Eff. 


dense core 


H2O maser 


22.2 


40 


50 


0.03 


Eff. 


outflow and/or disk ? 


CH3OH maser 


6.7 


130 


50 


0.1 


Eff. 


high-mass signpost 



Note. — Observatories: Eff.: Eff'elsbcrg, PV: Pico veleta, CSO: Caltech Submillimeter observatory 



Table 3 
Results 



source 


cm flux 


H 2 0. 


CH3OH 


MSX/cm MSX/mm 


cm/mm 


wings 


SiO 


CH3OH 


CH3CN 




[m.Iy] 


fJyl 


[Jy] 


["1 f'l 


["1 










05358 + 3543 


< 1 


45 


162. 1 






+ 


+ 


+ 




05490 + 2658 


< 1 


— 


- 


16.9 




+ 






- 


05553 + 1631 


1.3 


— 


- 


3.5 0.8 


2.8 


+ 


+ 


+ 


+ 


18089-1732 


0.9 


75 


54.0 2 


8.1 11.2 


3.2 


+ 


+ 


+ 


+ 


18090-1832 


? 


— 


82. 2 


2.6 




+ 








18102-1800 


44 




8.8 3 


1.7 25.3 


27.0 


4- 


4- 






18151 1208 




S* 


50. 4 






+ 


+ 


"f 




1 £i ]_ "j cj 1550 


in n 
' 


? 




6 Q 




+ 








18182—1433 


< l 


18 


24.0 


10.8 




+ 


+ 


+ 


+ 


18223— 1243 


< l 




~ 


3.6 




+ 








18247 — 1 147 


47 


~ 


\ 


7.4 5.4 


" 










18264— 1152 


< 1 


50 


3.8 1 


9.4 




+ 


+ 


+ 


+ 


18272— 1217 


110 


— 


— 


4.4 8.9 


5.3 










18290-0924 


7.0 


4 


15. 1 2 


12.7 4.5 


17.3 


+ 


- 


- 


- 


18306-0835 


82 


0.7* 


— 


15.5 11.6 


11.8 


+ 








18308-0841 


< 1 


1.5* 


— 


10.6 




— 


+ 


- 


+ 


18310-0825 


7.0 


— 


— 


5.1 16.8 


10.6 


+ 








18337-0743 


< 1 


— 


— 


4.7 




? 








18345-0641 


27 


3* 


5.4 4 


4.0 




4- 


4- 


4- 


4- 


18348-0616 


54 






14.8 18.6 


67.5 


+ 








18372-0541 


80 


1.5* 


9.2 1 


3.0 1.8 


1.9 


4- 








18385-0512 


29 


200 




4.4 2.8 


3.4 


+ 








18426-0204 


1.6 


7 


1.3 


6.7 12.0 


17.6 


+ 








18431-0312 


3.5 






30.3 2.4 


32.0 


+ 








18437-0216 


< 1 






6.9 




? 








18440-0148 


< 1 


? 


4.0 2 


2.6 




+ 


- 






18445-0222 


68 






0.4 18.5 


18.9 


? 








18447-0229 


1.6 






13.2 4.4 


10.0 


? 








18454-0136 


42 




2.5l 


4.6 


6.7 


+ 








18454-0158 


18 








2.2 


+ 








18460-0307 


< 1 






5.3 




+ 








18470-0044 


? 




7-5 3 


20.6 




+ 


+ 


+ 


+ 


18472-0022 


110 






2.1 


0.8 


+ 








18488 + 0000 


194 


1* 


16. 9 4 


19.0 4.1 


16.2 


+ 


+ 


+ 


+ 


18517+0437 


< 1 


45.3 


279 1 


22.6 




+ 








18521+0134 


< 1 




1.3 


3.7 












18530+0215 


311 






7.3 7.2 


3.8 


+ 


+ 


+ 


+ 


18540 + 0220 


97 






12.1 2.8 


10.3 


? 








18553+0414 


< 1 


50* 








+ 








18566+0408 


< 1 


3* 


7.2 3 


0.9 




+ 


+ 


+ 


+ 


19012+0536 


< 1 


2* 


0.9 


5.9 




+ 


+ 


+ 




19035+0641 


4.0 


9 


14. 3 5 


5.1 7.1 


9.0 


+ 


+ 


+ 




19074+0752 


14.8 






2.3 




+ 








19175+1357 


5.1 






2.1 


37.0 










19217+1651 


32 


9* 


0.9 




6.5 


+ 


+ 


+ 


+ 


19220 + 1432 


11.0 






4.8 6.1 


10.4 


+ 








19266 + 1745 


< 1 


1.5* 


l.Tl 


4.1 




+ 








19282+1814 


< 1 




2.9 1 


11.8 




+ 








19403 + 2258 


< 1 






23.5 




+ 








19410+2336 


1.0 


110 


30. 1 


3.8 




+ 


+ 


+ 


+ 


19411 + 2306 


< 1 






2.5 




+ 


+ 


+ 




19413 + 2332 


< 1 






11.1 




+ 


+ 


+ 




19471 + 2641 


? 










+ 








20051 + 3435 


< 1 






6.7 












20081 + 2720 


< 1 










+ 








20126+4104 


? 


15 


36. 6 


2.6 




+ 


+ 


+ 


+ 


20205 + 3948 


8.3 






12.1 11.5 


23.0 


+ 








20216+4107 


1.4 






4.3 




+ 








20293 + 3952 


7.6 


100* 




2.0 23.5 


15.0 


+ 


+ 


+ 


+ 


20319 + 3958 


25 






8.3 8.5 


15.5 


+ 








20332+4124 


136 


0.6* 






10.5 


+ 








20343+4129 


1.8 








13.2 


+ 








22134 + 5834 


3.7 








3.5 


+ 








22551 + 6221 


4.5 








74.7 










22570 + 5912 


29 






6.4 13.3 


7.3 


+ 








23033 + 5951 


1.7 


4 




6.4 10.6 


7.2 


+ 


+ 


+ 




23139 + 5939 


1.4 


400 


2.6* 


2.1 1.2 


1.0 


+ 


+ 


+ 


+ 


23151 + 5912 


< 1 










+ 








23545 + 6508 


1.0 










+ 








Note. — '* ': 


new dctcctio; 


is, '?': tc 


ntativc detections; the MSX/cm MSX/ 


mm and c 


m/mm colt 




present offsets between the different observations, and the last 4 column 


show detections 


( + ), non-dctcction (-) 


or non-observations (empty space); 1 


Szymczak et al. 2000 


; 2 Walsh ct al. 


1998; 3 Slysh ot al. 1999; 4 v. d. Walt ot al. 1995; 5 Caswell 



Vaile 1995; 6 MacLeod 1992 
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-10 -5 5 10 5 10 

[kpc] galactocentric distance [kpc] 

Fig. 1. — (a) The Galactic distribution of the sample: filled circles represent sources with no distance ambiguity, triangles show sources 
with distance ambiguities (filled triangles: near; open triangles: far). The size of the symbols scales with their luminosity. Squares show the 
Galactic center and the location of the sun (8.5 kpc from the center), while the circles indicate galactic distances in 2 kpc steps, (b) The 
Galactic scale height distribution of the sample. Symbols are the same as in (a). The horizontal lines show the average Galactic scale height 
of massive star formation regions 37 pc (Bronfman et al. 2000). 
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Fig. 2. — Examples of the HIRES images (contours) at 12 /irrt and 100 (im on grcyscale 1.2 mm dust images (Beuthcr ct al. 2001). While 
the bottom images are typical with both HIRES bands corresponding to the same mm core, the middle and top images are counterexamples 
with the 12 fim peaks splitting up or tracing different sources in the same field of view. 




Fig. 3. — Examples of HIRES 2 component greybody fits. While the dashed lines show the two greybodies separately, the full line shows 
the two component fit. Each panel presents on the top left the source name, luminosities and temperatures. 
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M gas [M ] M/L [M /L ] 

Fig. 4. — (a): Plotted are the derived FIR luminosities against the core masses. Triangles show HMPOs with distance ambiguity (filled: 
near distance; open: far distance) and filled circles without distance ambiguity. The crosses show data from UCHu regions (Hunter et al. 
1997, 2000). The grey shaded area defines the region in which 90% of the simulated clusters are found (Walsh ct al. 2001). (b): The 
histogram presents the distance independent quantity M/L with again a clear separation of HMPOs (grey) and UCHns (thick line). 
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Fig. 5. — Water (grey shaded) and methanol (thick lines) maser example spectra. The abscissa shows the velocities (LSR) for both spectr; 
while the ordinate shows the intensity of the water maser in [Jy], the methanol maser intensity is shown in arbitrary units. 



16 



2.5 



m 

o 



1 1 1 1 1 1 1 1 1 1 1 

a cm non— detections 






•* weak cm detections 






o methanol maser detections 


g 


- 


□ water maser detections 










- 


[®l m 












H a § 












© 






A- * • 






* A 






H 









0.5 



1 1.5 
log 10 [S(25) /S(12)] 



Fig. 6. — IRAS Color-color diagram of the whole sample: open triangles represent sources without and filled triangles with cm emission. 
The surrounding circles and squares represent detected methanol and water maser emission, respectively 
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Fig. 7. — The HIRES IR luminosities are plotted against the luminosities based on the cm observations. Again filled circles show sources 
with no distance ambiguities and triangles sources with distance ambiguities (filled: near distance; open: far distance) The straight line shows 
the 1:1 correlation. 



